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Methods 


Relaxation time constants can yield important biological information. Generating 
spatial parametric maps of T x and T 2 are well known MRI techniques. These maps 
usually involve fitting an exponential function to data points in the decay 
dimension for each image voxel. Chemical shift information is generally lost. 

In NMR spectroscopy, MRS and MRSI, it is known that various spectral peaks often 
have different relaxation times. Previous single-voxel-spectroscopy (SVS) studies 
of muscle in DMD have shown that the T 2 of the water signal changes in response 
to steroid treatment, and that the overall T 2 of the muscle changes with disease 
progression. 

Efforts to characterize the relaxation properties of magnetic resonance (MR) 
signals in the frequency domain have traditionally used either of two approaches: 

i) parametric line-shape modeling 

ii) spectral integration over the range of a peak of interest 


While pointwise relaxation fitting in the frequency 
domain is potentially broadly applicable, the present 
results are from data of limited scope. 

Acquisition: 

In vivo muscle data and phantom data acquired on a 
Siemens 3T whole-body MRI instrument with a 
STEAM SVS sequence. 

For vastus lateralis data, TE = (11ms, 27, 54, 243). 

TR = 9000ms. Half-echo FID sampling of 2048 points 
for 1024ms. Water placed on-resonance. 

Processing: 


These approaches involve possibly erroneous assumptions and may obscure 
useful information contained in the decay-dimension of a dataset. 

Here, we demonstrate fitting every point in the frequency domain for a transverse 
relaxation time constant (T 2 ) and equilibrium magnetization (M 0 ), from which we 
generate plots of T 2 and M 0 as a function of chemical shift. 


Custom python code incorporating apodization, peak 
alignment, automatic or manual phasing, zero-filling, 
and eventually pointwise relaxation fitting and 
plotting. Using existing libraries such as Numpy, Scipy, 
Matplotlib, pydicom, and NMR Glue. 
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Magnitude Real 





1 Hz exponential apodization on both magnitude and real 


Taking the magnitude of the complex-valued spectral data following FFT has some advantages: 

1) Phasing the spectrum to obtain a true "pure absorption mode" can be problematic, magnitude mode avoids this. 

2) The baseline can exhibit artifacts that cause relaxation fitting to fail or become unstable, magnitude mode mitigates this. 

3) In general, we have found highly stable results with magnitude mode fitting. 

However, magnitude mode may introduce spurious relaxation time fittings between peaks that arise from interference in the dispersive 
channel. 
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Spectral T 2 comparison between VL 
in a 14 year-old healthy male (top) 
and an 11 year-old male with DMD 
(bottom). 

The water T 2 at the center of the 
spectrum is quite similar between 
these subjects at approximately 
33ms, however the fatty-tissue 
replacement in the DMD subject 
reveals dramatic T 2 features at 2.7 
and 3.5ppm as well as T 2 variability 
among the lipid resonances 
between 0.7 and 2.5ppm. 
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VL spectral T 2 plots for a male 
with DMD, longitudinally 
tracking rapid disease 
progression. 

Only minor changes in the 
spectral T 2 are observed 
between 9 and 10 years-of- 
age. The prominent feature at 
approximately 3.5ppm is very 
clear and persistent in this 
subject. 
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Spectral T 2 plot from a 
vegetable oil phantom with 16 
TE values (11 - 288 ms). 

The cross marker overlays 
depict the T 2 values from 
subcutaneous adipose tissue 
at 7T as previously reported by 
Ren, et. al. 2 , and show 
remarkably close agreement 
for all lipid resonances with 
the exception of "C". The inset 
at lower left shows the voxel 
placement in the cylindrical oil 
and water (LEGO) phantom. 
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Spectral pointwise T 2 fitting was robust and reproducible, 
including at the tails of the spectral window where there 
is little obvious signal. The T 2 of the dominant spectral 
component, often ^O, determined the mean T 2 values at 
the edges of the spectrum and indicate the broad extent 
of signal contribution. Variance of the T 2 fit increases with 
spectral offset from the signal maximas. 

VL 1 H 2 0 (4.7ppm) T 2 values were found to occur in the 
narrow range of 26 - 39ms, with DMD afflicted muscle 
showing a greater range that tended to decrease with 
disease progression (Figure 3). The lipid methyl resonance 
at 0.9 ppm had a T 2 range of 55 - 90ms, and generally 
lower T 2 than the methylene signals at 1.4 and 1.7 ppm. 
Creatine and trimethylamines are not well resolved but 
may contribute to the spectral T 2 peaks at 3.0ppm and 
3.2ppm. In the higher fat-fraction VL spectra, a large and 
broad T 2 feature appears at 3.5ppm to 4.0ppm and 
warrants future investigation. 

Analysis of the oil and water phantom spectra shows 
highly similar T 2 values for 6 lipid resonances previously 
studied (Figure 4). 2 


The ability to visualize a T 2 spectrum aids in interpretation of complex signal 
behavior, and identification of overlapping signals. One motivation for 
pointwise T 2 fitting was to better understand the spectral heterogeneity of 
1 H 2 0 T 2 in DMD, especially later in the disease when lipid is the dominant 
signal. At high lipid fractions, the apparent water T 2 in muscle often 
decreases, and the contribution from the proximal lipid resonances does 
not appear to be the cause. The spectral T 2 provides a unique fingerprint 
that conveys additional information on muscle pathology. 

Future directions include a rigorous assessment of the relative benefits of 
fitting on real or magnitude spectra. The TE sampling scheme will be 
investigated as it may play a role in some of the spectral T 2 features. 
Additionally, inversion recovery SVS acquisition can be used to generate a 
spectral T v A bi-exponential decay model may yield additional information 
and help to distinguish otherwise hidden spin populations. More work is 
required to link specific T 2 features to normal and pathological biology. 

These derived spectral relaxation plots may lead to unprecedented richness 
in relaxation analysis. This approach can readily be applied to other tissues 
in vivo, NMR spectroscopy in vitro, and coherent spectroscopy in general. 
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